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Abstract 

We study the available data in polarized e — p deep inelastic scattering to test two different 
solutions to the so called spin crisis: one of them based on the axial gluon anomaly and 
consistent with the Bjorken sum rule and another one, where the defects in the spin sum 
rules and in the Gottfried sum rule are related. In this case a defect is also expected for 
the Bjorken sum rule. Experimental data, especially the very recent SLAG E154, favour the 
first solution and demand a gluon polarization AG = 2.25 it 1.39. 
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1 Introduction 



The earlier EMC CERN experiment |j^] and the importance of testing the Bjorken sum rule 
[^] have stimulated a considerable experimental activity in measuring g^, {gi^^) and gf. 
The result found from EMC for the first moment, 

Tl= I = 0.126 ±0.010 ±0.015, (1) 

JO 

has been confirmed by the SMC CERN experiment Q (at <Q^>= 10 GeV"^ , which is almost 
the same as for EMC) and SLAC 0] (at <Q'^>= SGeF^), giving respectively 

V\ = 0.136 ±0.011 ±0.011 {SMC), 
V\ = 0.127 ±0.004 ±0.010 (-E143). 



(2) 



These experiments have also measured the deuteron structure function |^, ^ and from Pf 
by substr acting Fj", one gets 



(3) 



T^l = -0.063 ±0.024 ±0.013 {SMC) 
= -0.037 ±0.008 ±0.011 (S143). 

At SLAC with polarized He^ targets they also obtained Q 

= -0.031 ± 0.006 ± 0.009 (E142, = 2GeV'^), (4) 

and the preliminary result 

/ gi{x)dx = -0.037 ± 0.004 ± 0.010 (-E154, = ^CeV"^). (5) 

J0.014 

The preliminary result from Hermes [| (T^ = -0.032 ± 0.013 ± 0.017 at <Q'^>= ^CeV'^) 
is consistent with SLAC data. 

The main issue of this experimental work is to test the validity of the Bjorken sum rule, 
which up to O(as^), for nj = 3, is given by ||l^ 

2 / \ 3" 



(6) 



Indeed, an interpretation of the defect in the Ellis and Jaffe sum rule |11] for implied 
by Eqs. (|I|) and (§) has been given in terms of a negative contribution coming from a large 
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positive polarization of the gluons AG [12|, which is the same for proton and neutron, such 
that it does not affect the Bjorken sum rule. 

An analysis of the existing data, excluding the most recent and very precise SLAC E154 
data, has been performed in a framework consistent with the Bjorken sum rule and including 
next to leading order (NLO) effects in the evolution equations to relate data at different Q^, 
and it provides a fair description of the experimental results |13|. 

All the existing data do not exhibit a clear evidence of evolution, i.e. within error 
bars they are compatible with scaling behaviour in the entire range accessible by all 
experiments. Although a complete NLO analysis of the g\ data is the correct procedure 
without a doubt, we are aiming to demonstrate that the most accurate data provided by 
the SLAC experiment is really telling us something important about the gluon contribution, 
even at a level of a less sophisticated leading order (LO) analysis. In spite of the fact that 
the SMC CERN experiment can reach higher and smaller x, the actual precision achieved 
cannot provide a reliable test for the NLO theoretical analysis. 

Here we want to compare the current interpretation of the defect in the Ellis and Jaffe 
sum rule in terms of a large flavour singlet contribution to the nucleon polarization coming 



from the gluons, with another one [14|, where one relates this defect to the one in the 
Gottfried sum rule [15| and to the role that Pauli principle seems to play |16]. This is done 
by relating the first moments and the shapes of the parton distributions, as first proposed 
in Ref. [0]. Indeed, if one assumes the validity of the Adler sum rule |17|, 

u — d = [(n — n) — (d — d)] + -u — d = 1 + u — d, (7) 

the defect in the Gottfried sum rule implies 

u-d<\. (8) 

If one thinks that the Pauli principle is responsible for the inequality (P), it is reasonable 
to assume that it is u^, the most abundant valence parton, which receives less contribution 



from the sea, so that we have |14] 

Alt = v} — v}- ^Uyal +u — d, (9) 

producing a defect in the Ellis and Jaffe sum rule for the proton 

Ar? = ^ (n - d) = ^ (-0.15 ± 0.04) ~ -0.033 ± 0.009, (10) 
9 9 
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in fair agreement with the experiments. 

An empirical test for the two interpretations might be given using the experimental 
information on the x-dependence of the polarized structure functions. The higher preci- 
sion of SLAC data (especially the ones of E154 for g") and the agreement between E142 
(<(5^>= 2GeV'^) and E154 {<Q'^>= BGeV'^) data, suggest to describe them together with 
the E143 measurements {<Q'^>= 3GeV^) on proton and deuteron in terms of the same 
parton distributions, and we consider for these distributions two options corresponding to 
the two different solutions to the spin crisis. 

We shall neglect higher-twist terms, supported by more recent theoretical evaluations of 



the contribution of these terms [ 18 1 which lead to results smaller in modulus and sometimes 
opposite in sign than the previous one, consistent with an experimental determination of 
these terms by the SLAC group [^] . 

The paper is organized as follows. In the forthcoming section we shall describe the SLAC 
data, with proton and deuteron targets at <Q^>= 3 GeV^ and He^ target at <Q^>= 2 GeV^ 
(E142) and at <Q^>= 5 GeV'^ (E154), with the two different options. Then we shall present 
the method we used to solve the Altarelli-Parisi evolution equations and to find the parton 
distributions at = 10 GeV"^ which we shall compare with CERN data. Finally, we shall 
give our conclusions. 

2 Description of SLAC data 

We describe the proton and neutron polarized structure functions at Qq = 3GeV^, in terms 
of the valence quark and gluon polarized distributions only, using a simplified version of the 
functional forms used in Ref.pO|] (in our case we take 7g(g = u, d, G) =0), namely 

xAuy{x,Ql) = rjuAux"^ {1 - xf" , 

xAd^ix,Ql) = ris,AdX^^(\-xf^, (11) 
xAG(x,Ql) = r]GAGx''o^i_^)ba^ 

where rjq {q = u, d, G) are the first moments of the distributions and Aq = Aq{aq, bq), 

+ (12) 

^ Jo T[aq + bq + l) 

in such a way that 

/ dxAqx''''-\l - xp = 1. (13) 
Jo 
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As pointed out by several authors |21|, to avoid the inclusion of soft contributions into the 
coefficient functions one has to choose a factorization scheme in which the gluon polarization 
contributes to the first moments of and 5" (for Uf = 3) : 



(14) 



The gluonic term appears to be a higher order correction but is not, because tig{Q'^) rises 
logarithmically with and, if the gluons had a positive polarization, it could, in principle, 
be large enough to explain the defect in the Ellis and Jaffe sum rule. 
We take 



af'^Hx.Q') = ^(3^)An.(x,Q2) + lQ)A4(x,Q^) 



6 TT 



(Act ® AG) (x, 



gfix,Q' 



(15) 



where ujd = 0.058 |2^] takes into account the small D-wave component in the deuteron 
ground state. In Eq. ( [l5| ) the QCD corrections in the quark sector are included in the F 
and D values entering in the expressions of the first moments of the quark distributions 
(F = 0.46 ± 0.01, D = 0.79 ± 0.01 ||): 

2 



Dm 



5F 1 



TT 



(10.46F + 2.48L>) 



IT 



20.22 (2F + D) 



as(Q') 



TT 



5D 1 



TT 



20.22 {3F + AD) 



(7.44F + 15.421?) 



(16) 



vr 



and the gluon contribution appears as a convolution |24] 

(Act® AG) {x,Q^) = 
We fix (as (3 GeV"^) = 0.35 ± 0.05) 



— (1 — 2 zj In 



1) ag(^,q2 



UQl) = FiQl) - D{Ql) 



-0.26 ±0.02, 



(17) 



(18) 
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and we explore the two options A and B, the first one with 

iquiQl) = '^F{Ql) = 0-76 ± 0.04, (19) 

and riG free, the second one with r/^ free and r]G = 0. Options A and B correspond respec- 
tively to the interpretation of the defect in the Ellis and Jaffe sum rule for T\ in terms of 
the anomaly, assuming that the Bjorken sum rule is obeyed, and to the case of a smaller Au 
resulting from the Pauli principle. 

Since we know that v) dominates at high x and that the gluons dominate in the small x 



region, we restrict, as in Ref.||25[], the values of the parameters in Eqs. (|llD, to be consistent 



with the information we already have for the parton distributions, by the following limitations 

6„>1, 6rf>3, 6g>5, (20) 

and we assume 

ttu = ad- (21) 
Indeed, especially for option A, where one describes two functions, g^{x) and gi{x), in 



terms of three distributions, Eqs. (11), one has to make sure to exclude some choices of the 
parameters describing well the data, but not consistent with the information one has from 
the unpolarized data, that is, e.g. about 1/2 of the proton momentum (in the = oo 
frame) is carried by the gluons and that the partons are dominating the high x region. 

The parameters corresponding to the best fit of the SLAG proton and deuteron data 
for options A and B are given in Table while in Figs. 1, 2 and 3 one compares the two 
resulting curves with SLAG data. 

All the data are well described with the two options, except for the ones by E154, which 
are better described by the option A, the one with gluon contribution and consistent with 
the Bjorken sum rule. Option A implies a large value of AG = 2.25 it 1.39. 

3 Parton evolution equations 

For the polarized parton distributions one has the Dokshitzer-Gribov-Lipatov-Altarelli-Parisi 
evolution equations (DGLAP) |^6|, which are, in the variable t = IuQ^/Aq^^, and at LO in 
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Us (Ag* = xA(7* and = xAG), 



dt 



-{x,t) 



dAg 
dt 



+ 



xA) 



2tt [A 3 



dz (z-\]A~g(-,t 



(i = l,...,2n/) 



27r 



1 /I ^"-^ 
^ dz-{2-z)Y,Ait(^^.t 



+ dz3 



A5(-,t). 



(22) 



We work in the fixed flavour scheme where the number of flavours in the sphtting functions 
is = 3, while for the Q^-evolution of Ug one has 

127r 



(33 -2n/) In-S 



2 ' 



(23) 



A, 



(4) 2 
QCD 



with aJ^^^ = 201 MeV to reproduce, with n/ = 4, as{3GeV'^) = 0.35 ± 0.05. 



For the solution of the DGLAP equations we have used a method |27] that consists in 
expanding the parton distributions into a truncated series of Chebyshev polynomials, 

r) 1 n— 1 

p^(x,t) ^p^(T(x),t) = - ^ ^ pX^s,t) r^r,) r,(r(x)), (24) 

s=0 1=0 

where Ti are the Chebyshev polynomials and 



t{x) 

Umax 
VI 



2 In X Umax 
Umax 

4 In 10, 

f 0.5 Z = 
I 1 Z>0 ' 



(25) 



exp 



(Ts + 1) 



2 s + 1 

cos vr. 

2 n 



Substituting Eq. 

dpi 
dt 



in Eqs. ( |2^ ) gives rise to a system of coupled differential equations 



%t) 



2tt 



7 n-1 



E J2^ksl^sit), (^ = l,-,7; fc = 0,...,n-l). 



(26) 



j=l s=0 
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in which = p^{xk,t) are the values of the polarized distributions Ag* and Ag in the n 
points Xk corresponding to the nodes of the Chebyshev polynomials. 

With the initial conditions given by the results of the fits A and B to the SLAC data we 
get the evolved distributions at = WGeV^. 

4 Comparison of the evolved distributions with experiments 

The predictions for the evolved distributions at = WGeV^ are compared with CERN 
measurements j^, ^ at = 10 GeV^ in Figs. 4 and 5. There is a better agreement for 
option B for the proton (total of 10.6 for 12 points to be compared with = 16.2 for 
option A) while for the deuteron the option A has a slightly lower (total of 14.2 for 
12 points to be compared with x^ = 16.8 for option B). Note that for option A one has six 
free parameters, but only four for option B. 

It is interesting to remark that with both options one fails to reproduce the rise at low 
X of xgi{x), which turns negative below x = 10^^ for option A, while xgf{x) in the same 
region is in agreement with the trend of the data (see Fig. 5). This is due to the isoscalar 
nature of the anomaly contribution, which is expected to be the same for proton, neutron 
and deuteron (neglecting the small correction coming from the D-wave component in its 
ground state). 

5 Conclusions 

We have studied the precise SLAC data on polarized structure functions with the main 
purpose of testing the Bjorken sum rule and the necessity of a gluon contribution to the 
polarized structure functions. As a result we find a better description of the data, especially 
of the very precise ones by the SLAC E154, with a gluon contribution and imposing the 
Bjorken sum rule than for the option without gluons and with the first moment of Au free. 
The best fit in the first case is obtained with a rather large value of AG = 2.25 it 1.39. 

The AG distribution appears to be more singular than the Ag's, leading to the conclusion 
that gi{x) should also become negative at small x. These results are in agreement with what 
was found in Ref. [^]. Although we have found here ac = 0.44 and ac = 0.13 in [^], these 
two different powers lead to comparable values for /q q\ AG(x) dx, which is 1.5 in Ref. p^ ] 
and 1.0 here. Note that in Ref. [^^, since the (1 — x) terms were omitted, one expects to 
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get, for gluons and for valence quarks, a slightly more singular behaviour when x ^ than 
in the present analysis. 

Concerning the test of the Bjorken sum rule, the fact that one gets a better fit by allow- 
ing the presence of a contribution coming from the gluons speaks in favour of its validity. 
Assuming the validity of the Bjorken sum rule, this contribution was advocated to explain 
the defect of the Ellis and Jaffe sum rule for the proton. However, the fact that, according 
to our description, at low x the isovector contribution is expected to be overwhelmed by the 
isoscalar one, which is more singular, suggests that a precise test of the Bjorken sum rule 
rather requires very accurate measurements in the x region where the two contributions are 
comparable. Indeed, in the very small x region, where gi{x) — gi{x) should have a small 
value being the difference between two almost equal negative large quantities, normalizations 
uncertainties can produce substantial errors. SLAC E154 has measured gi{x) with an out- 
standing precision and we look forward to have, in the same x range, a comparable precision 



for gl from SLAC E155 |29|. 
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TABLE I 





A 


13 
JD 


au = ad 


0.79 ± 0.05*^ 


1.03 ± 0.07*^ 


bu 


1.51 ± 0.17'- 


1.88 ± 0.23'- 


bd 


3.00 ± 0.13(*) 


4.52 ± 0.52^*) 


aa 


0.44 ± O-IS*- ^ 




ba 


5.00± 1.05(*) 




Vu 


0.76 ±0.04 


0.63 ±0.03^*) 


m 


-0.26 ±0.02 


-0.26 ±0.02 


VG 


2.25± 1.39(*) 




Xndf 


0.84 


1.09 



Table 1: The results of the options A and B for the values of the parameters of the fits at 
= 3GeV^. The free parameters are marked with an asterisk. 
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Figure Captions 

Fig. 1 The best fit for the options A (sohd line) and B (dashed line) (see text) are compared 
with the SLAC data on proton for xg^ix) at <Q'^>= 3GeV'^ from Ref. [|]. 

Fig. 2 Same as Fig. 1 for the deuteron SLAC data for xgi{x) at <Q^>= 3GeV^ from 
Ref. §. 

Fig. 3 Same as Fig. 1 for the neutron SLAC data for xgf{x) at <Q'^>= 2 GeV"^ from Ref. 
(triangle) and at 5GeV^ from Ref. [^] (boxes). 

Fig. 4 The data on proton for xg[{x) from SMC at <Q'^>= lOGeV"^ from Ref. are 
compared with the results of the options A (solid line) and B (dashed line), evolved 
to Q2 = lOGeV'^. 

Fig. 5 Same as Fig. 4 for the deuteron SMC data for xgf{x) from Ref. Q. 
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